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We analyze general properties of the perturbation expansion for two-dimensional quantum critical 
metals with singular forward scattering, such as metals at an Ising nematic quantum critical point 
and metals coupled to a U{1) gauge field. We derive asymptotic properties of fermion loops 
appearing as subdiagrams of the contributing Feynman diagrams - for large and small momenta. 
Substantial cancellations are found in important scaling limits, which reduce the degree of divergence 
of Feynman diagrams with boson legs. Implementing these cancellations we obtain improved power¬ 
counting estimates that yield the true degree of divergence. In particular, we find that perturbative 
contributions to the boson self-energy are generally ultraviolet convergent for a dynamical critical 
exponent z < 3, and divergent beyond three-loop order for 2 : > 3. 
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I. INTRODUCTION 

The low-energy behavior of interacting fermion systems 
can be strongly influenced by collective bosonic degrees of 
freedom such as critical order parameter fluctuations or 
emergent gauge fields. In metallic systems, the scattering 
of electrons by gapless bosons destroys Landau quasi¬ 
particles and thus leads to a breakdown of Fermi liquid 
theory [1]. These effects are particularly pronounced in 
two-dimensional systems, and they have therefore fre¬ 
quently been discussed in the context of layered com¬ 
pounds such as cuprate- or iron-based high-temperature 
superconductors. In this paper we focus on systems in 
which the critical bosons carry a small momentum and cou¬ 
ple to the electron charge, leading thus to a singular but 
spin-conserving forward scattering of electrons. This class 
of systems comprises electron liquids coupled to emergent 
U{1) gauge helds [2], as well as quantum critical metals at 
the onset of nematic or other translation-invariant charge 
order [3-5]. 

The breakdown of Fermi liquid behavior due to singular 
forward scattering is revealed already by first-order pertur¬ 
bation theory. One-loop results for the boson and fermion 
self-energies in two-dimensional metals were first derived 
for fermions coupled to a 17(1) gauge field [6], and later for 
systems at a nematic quantum critical point (QCP) [3,4]. 
The bosons receive a Landau damping term generated 
by fermionic particle-hole excitations, and the one-loop 
fermion self-energy scales as [w]^/^ at low frequencies w. 
This self-energy dominates over the linear frequency term 
in the fermion propagator and thus implies a destruction 
of Landau quasi-particles. The main contributions to 
the fermion self-energy at a fixed Fermi momentum 
come from particle-hole excitations near the same kp 
and its antipode — k^?, with small momentum transfers q 
tangential to the Fermi surface and excitations energies 
of the order jqj^. The coupled fermion-boson theory is 
scale-invariant with a dynamical exponent z = 3 at the 
one-loop level [7]. 

An early analysis of two-loop corrections did not re¬ 


veal any qualitative changes of the boson and fermion 
propagators [8]. The loop expansion was expected to be 
controlled by the inverse fermion flavor number iVy, as 
usual [7]. However, the whole case was opened again when 
S.-S. Lee [9] discovered that the naive 1/Nf expansion 
breaks down and Feynman diagrams of any loop order 
contribute even in the limit Nf —?> 00 . Subsequently, 
Metlitski and Sachdev [5] constructed a general scaling 
theory for two-dimensional metals at a nematic QCP and 
the formally similar C/(l) gauge field problem. Symme¬ 
try constraints restrict the theory such that only two 
independent anomalous scaling exponents are possible: 
an anomalous dimension of the fermion fields rjf and an 
anomalous dynamical scaling exponent z 3. A small 
contribution to rjf was indeed found at three-loop order, 
while the dynamical exponent remained unrenormalized 
at z = 3 at that loop level [5] . Most recently it was shown 
that the absence of renormalizations of z up to three-loop 
order has special reasons that do not apply to all higher- 
order contributions. A divergence leading to anomalous 
dynamical scaling was found in a class of four-loop dia¬ 
grams [10]. Surprisingly, the divergence turned out to 
be stronger than logarithmic, indicating an anomaly or 
instability whose nature has not yet been clarihed. 

Naive power-counting frequently overestimates the ac¬ 
tual degree of divergence of perturbative contributions. 
Cancellations may occur for single Feynman diagrams, 
due to oscillating integrands, and also between distinct 
Feynman diagrams. In this context, the scaling behav¬ 
ior of fermion loops plays a particularly important role. 
Most Feynman diagrams contain fermion loops that are 
connected to each other and/or to open fermion lines 
by bosons. All hitherto observed cancellations are re¬ 
lated to cancellations within a loop or between loops with 
permuted vertices. 

In this paper we analyze the asymptotic behavior of 
fermion loops (with N vertices) and sums of loops with per¬ 
muted vertices in the scaling limit that applies to fermions 
coupled to bosons in two dimensions as described above. 
We perform the analysis for a generalized class of theories 
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with dynamical exponents z > 2 , as originally introduced 
by Nayak and Wilczek [11]. Systematic cancellations are 
obtained if one or several of the external loop momenta 
are much smaller than the others. As an application, 
we can assess the degree of divergence of several classes 
of Feynman diagrams without embarking on a tedious 
specific calculation of the corresponding integrals. 

The paper is structured as follows. In Sec. II we intro¬ 
duce the low-energy effective quantum field theory for the 
systems described above, and we provide the precise defi¬ 
nition of the fermion loops. The result of the elementary 
integration over the momentum in the loop is presented 
in Sec. III. The asymptotic behavior of fermion loops 
with three vertices is analyzed in Sec. IV. An important 
application is a deeper understanding of the mechanism 
leading to the absence of a renormalization of the dynam¬ 
ical exponent z in the three-loop calculation by Metlitski 
and Sachdev [5]. Asymptotic properties of fermion loops 
with an arbitrary number of vertices N are derived in 
Sec. V. In particular, it will be clarified to what extent 
cancellations obtained at three-loop order can be expected 
also at higher orders. In Sec. VI we provide several exam¬ 
ples of improved power-counting. We finally summarize 
the main results in Sec. VII. 


II. FIELD THEORY AND FERMION LOOPS 


self-energies, the boson and fermion propagators are ob¬ 
tained as [ 6 ] 

= ( 2 ) 

Gs^W = skx + ky - (3) 

where q = (go, q) and k = (fco, k) contain momentum and 
frequency variables, and 7 and k are positive constants. 
These propagators solve the RPA equations also self- 
consistently [7]. The linear frequency term in the fermion 
propagator is subleading compared to the fermion self¬ 
energy and has therefore been discarded. Note that D{q) 
does not depend on q^. 

In the following we consider a generalization of the the¬ 
ory where qy is replaced by qy~^ in the boson propagator, 
such that 

= ( 4 ) 

This generalization was introduced by Nayak and 
Wilczek [11] for the sake of an expansion in e = z — 2. 
More recently, it was used to define a manageable large- 
Nf limit of the theory [12]. For any z > 2, the one-loop 
fermion propagator then assumes the generalized form 


The low-energy behavior of two-dimensional fermion 
systems coupled to critical bosons with small momenta can 
be described by an effective field theory whose fermionic 
states are restricted to two Fermi “patches” near fixed 
Fermi momenta kp- and —kp- [5,9]. Choosing momentum 
variables such that the a:-component is perpendicular and 
the y-component tangential to the Fermi surface at ±ki?, 
the Lagrangian of the field theory can be written in the 
form [5] 


L = ^ - isdx - dy) 'ips 

S —± 

- 9s4’'^li^S - ^{dy4>f. ( 1 ) 

s=± 

Here 0 is a bosonic scalar field, while '0±, ili\. are fermionic 
fields with Nf flavor components corresponding to states 
with momenta on the the two patches near Tk^?. In the 
U(l)-gauge field problem, </> is the transverse gauge field 
and g+ = —g-. For the Ising nematic QCP, </> is the 
order parameter field and g+ = g-. In both cases the 
physical flavor number is A^/ =2. The derivatives are with 
respect to real space and imaginary time variables. Several 
numerical prefactors have been absorbed by a rescaling of 
fields and space coordinates. In particular |ys| = I. The 
gradient term in x-direction {dx4>)'^ is irrelevant under the 
scaling that emerges on one-loop level, and has therefore 
been discarded in the Lagrangian. 

In random-phase approximation (RPA), which corre¬ 
sponds to a one-loop calculation of the boson and fermion 


Gs{k) 


1 

skx + fcy - i{ko} ’ 


( 5 ) 


where 


The fermion propagator scales homogeneously as 
G{k) I—>■ X~^G{k) under the scaling of momentum and 
frequency variables 

ky I —y Xky^ kx '—^ X‘^kx^ kf) 1 —y (7) 

The boson propagator scales as D{q) >-)• X^~^D{q) under 
the analogous scaling of yoi <lx, and qy. For z = 3 we 
recover a = 1/3, and Digq) scales with the same power 
as G{k). 

Higher-order contributions (beyond one-loop) to the 
fermion and boson self-energies are analyzed in an expan¬ 
sion around the one-loop fixed point, that is, by comput¬ 
ing Feynman diagrams with one-loop propagators [5,8,9]. 
Note that this procedure can be viewed as a simple shift 
of the expansion point, where the one-loop self-energy is 
added to the “bare” part of the action and subtracted as 
a counterterm from the interaction part. 

The computation and asymptotic analysis of higher- 
order contributions is done most efficiently by evaluating 
first the fermion loops in the Feynman diagrams, and then 
the remaining bosonic momentum integrations. The N- 
point fermion loop on patch s is defined as the integrated 
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Figure 1. Graphical representation of the A^-point fermion 
loop. 


product of N fermion propagators, 


nAr,s(gi, ■ • •, Qn) = In,s{pi, ■ ■ ■ ,Pn) 





N 

l[Gs{k-p,). 

2=1 


( 8 ) 


The variables pi and qi are related by 


Qi = pi+i -p^ for i = -I 

Qn = Pi - Pn ■ (9) 

Note that I^^siPi, ■ ■ ■ ,Pn) is invariant under a global 
shift of momenta Pi Pi + q, and qi + ■ ■ • + qN = 0 due 
to energy and momentum conservation. Il^^siq!, ■ ■ •, 9Ar) 
can be represented graphically as a fermion loop, as shown 
in Fig. 1. 

The sum of Feynman diagrams contributing at a certain 
order in the loop expansion can be written in terms of 
symmetrized fermion loops. 


V 


with a sum over all permutations oi qi,... ,q]\[. The pref¬ 
actor 1/A^ compensates the TV-fold multiplicity arising 
from cyclic permutations. In various limits the sym¬ 
metrized loops are much smaller than the unsymmetrized 
contributions, due to systematic cancellations. Strong 
cancellations between contributions to symmetrized loops 
are well-known in one-dimensional Luttinger liquids [13], 
and have also been established for two-dimensional Fermi 
liquids [14,15]. 

In the remainder of our paper we evaluate the A^-point 
loops and analyze their asymptotic behavior in various 
important limits relevant to quantum critical metals, from 
which we can derive improved power-counting estimates of 
Feynman diagrams. To simplify the notation, we specify 
to fermion loops defined on the Fermi patch with s = -I- 
and drop the patch index. Results on the other patch 
(s = —) can be obtained from those for s = -|- by simply 


switching the sign of the cc-component of all momenta. 
Furthermore, we absorb the prefactor n/Nf in the fermion 
propagator by a rescaling of the frequency variable. Hence, 
{fco} stands for fco/|A:o|“ in the following sections. 


III. INTEGRATION OF FERMION LOOPS 


The fermion loop as defined in Eq. ( 8 ) involves a mo¬ 
mentum and a frequency integration. The momentum 
integral can be performed analytically by using the residue 
theorem (see Appendix A), yielding 


In{pi, 




i<j 


>PiO 


^Piy ~ Pjy ' 


X (jpiy Pjy) ^ Jijliko) (11) 


for any N > 3. The indices i,ji I mn from 1 to A with 
restrictions as indicated, 0 is the step function, and 


with 

k^ijl Pix(,Ply Pjy) ~b ; (^8) 

Fiji = {Pjy ~ Piy){piy ~ Pjy){Piy ~ Ply) > (14) 

^^Jl{ko) = {fco -Pio}{piy -Pjy) + cyc, (15) 


where “eye” denotes cyclic permutations of the indices 
i,jjk Note that the remaining frequency integration is 
limited to finite intervals and hence convergent. The 
fermion loop integral does not require any ultraviolet 
regularization. In the static limit qio = 0, all frequency 
variables pio are equal, such that In{pi, ■ ■ ■ ,Pn) = 0. 

The quantities Diji, Fiji, and ^liji{kf)) are antisymmet¬ 
ric in the indices i,j, 1. Hence, Jijiikg) is also antisymmet¬ 
ric in its indices. Using this antisymmetry, the three-point 
loop can be written in the particularly simple form 


h{Pl,P2,P3) 


E 

(ij)=(1.2).(2.3).(3,l) 


r'° dko q/ Pm -P 30 \ 
jp,0 47r yp^y-Pjy' 


X Jl23(fco) • 


(16) 


Under the scaling Eq. (7), the quantities Diji, Fiji, and 
Uiji(fco) scale homogeneously as 


Fiji I— > }?Diji, Fiji ^Fiji, ^iji^ko) I— > X^niji{ko). 

(17) 

Hence Jiji^ko) i—>■ A ^Jiji^ko) and 

In{pi, - ■ ■ ,Pn) In{pi, ■ ■ ■ ,Pn) , (18) 


in agreement with direct power-counting applied to the 
definition oi In{ pi,... ,pn) as a loop-integral. The sym¬ 
metrized fermion loops n)^“(gi,..., g^r) exhibit the same 
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scaling behavior. As a consequence, effective boson inter¬ 
actions with amplitudes proportional to 11 ^™ (< 71 ,..., 
are marginal for all N [16]. 

In the special case a = 0, corresponding to a bare 
fermion progagator, the fco-dependence in Vliji{kQ) cancels 
in the cyclic sum, such that — fliji = pio{pjy — 

Ply) + eye. The frequency integration in Eq. (11) then 
becomes trivial, yielding 

I%{pi,...,Pn) = — —){P^o-Pjo) 

X {piy — Pjy) ^ Jijl ■ (19) 


In particular, the three-point loop with bare propagators 
is obtained as 


l3iPl^P2,P3) 


1 

47r 


E 

(i,j) = (1.2).(2.3).(3,l) 


/ Pz0-Pj0 \ 
^Piy — Pjy ' 


X {Pio — Pjo) J123 ■ (20) 

In the infrared scaling limit pio 1 —>■ X^pio, pix X^Pix, 
Piy I —Xpiy with z > 2 and A — >■ 0, the term IIi 23 is 
subleading compared to D 123 and ^123 in the denominator 
of Ji 23 - In that limit, Eq. (20) is consistent with the result 
for the 3-point loop with a full Fermi surface integral 
obtained previously in Ref. 16. 


IV. ASYMPTOTIC BEHAVIOR OF 
THREE-POINT LOOP 


In this section we derive asymptotic scaling properties 
of the three-point loop. In particular, we analyze the 
behavior in a limit that determines the ultraviolet behavior 
of the Aslamasov-Larkin-type diagrams depicted in Fig. 2. 
These diagrams have been computed explicitly in the 
static limit (external gg = 0) by Metlitski and Sachdev [5] . 

Denoting the external frequency-momentum variable 
in the Aslamasov-Larkin diagrams by g, and the bosonic 
loop integration variable by g', the contributing fermion 
loops are n 3 (g, g', -g-g') = 73 ( 0 , g, g+g') and n 3 (g, -g- 
g', g') = 13 ( 0 , g, -q') = 13 ( 0 , -g, -g - g'). In the static 
limit go = 0, Eq. (16) yields 


/o 47r 
> 4123 (^ 0 ) j 


*(0,e(^)-e(^) 

7o 47r [ \qy + qU \qU 


( 21 ) 



Figure 2. Aslamasov-Larkin type diagrams contributing to 
the bosonic self-energy in a coupled theory involving fermions 
(solid lines) and bosons (wiggly lines). 


with 

Jl23{ko) = [q'xdy - QxQy - dydyidy + Qy) 

+i {{ko} - {ko - q'g}) qy] \ (22) 

Hence, 13 ( 0 , g, g-hg') with gg = 0 is non-zero only if jg^| < 
jgyl and sgn(g(^) = —sgn(gy). This kinematic constraint 
restricts the q'y integral to an interval of length jg^j. 

Now we can assess the ultraviolet (UV) behavior of 
the g'-integral in the Aslamasov-Larkin diagrams. We 
introduce an ultraviolet cutoff A such that [gg] < A^, 
jg^l < A^, and jg^j < A. However, for fixed g^, the q'y- 
integral is effectively restricted by g^. With this restriction, 
one has 


Di23 ^ qyX)^, Fi23 ^ Qy, A2l23(fco) Qy^^ , (23) 

and thus Ji 23 {ko) ^ (g^A^)”^, such that the UV behavior 
of the three-point loop is 

43(0,g,g + g') ~ -■ (24) 

qy 

The Aslamasov-Larkin diagrams contain two three-point 
fermion loops and two boson propagators. For g(^ ^ qy 
and gg ^ A^ the latter scale as qy/h.^. The integration 
measure scales as A^A'^qy = A^+^g^. Hence, the integral 
diverges in the UV limit as gj,A*“^. This agrees with 
the result obtained in a more explicit calculation for z = 
3 by Metlitski and Sachdev [5]. Note that the above 
contribution is obtained only if the two fermion loops are 
integrated over antipodal Fermi patches. For loops on 
the same Fermi patch, the g(,-integral vanishes since the 
integration contour can be closed without encircling any 
poles in that case [5,9]. 

The sum of all Aslamasov-Larkin diagrams can be ex¬ 
pressed by symmetrized loops. The symmetrized three- 
point loop can be written as 


n3 ^"'( 9T 72 , 93 ) = l 3 iPl,P 2 ,P 3 ) + I 3 {-Pl,-P 2 ,-P 3 ) ■ 

(25) 

Under the inversion pi 1 — > —pi and fcg 1 —> — fcg, the quanti¬ 
ties Diji, Fiji and Hij;(A:g) transform as 


Fiji ' t Diji^ Fiji ' t Fijl^ ^ijli^ko) I ^ ^ijl)k()) . 

(26) 

The symmetrized three-point loop contributing to the 
Aslamasov-Larkin diagrams in the static limit can thus 
be written as 


/f”(0,g,g + g') 


qo=o n° dko 

0( 

7^ ^ 

1 

© 

■ 0 

Jo 47r 

\ 

Wy + 7y^ 



1 

Di23+Fl23 + iFi23{ko) 


1 

Di23 — 7^123 +*I2i23 (^o) 


27) 


In the ultraviolet limit the leading terms in the last set 
of brackets cancel, leaving a much smaller contribution of 
order g^/A^, such that the symmetrized three-point loop 
scales as 


43^“(0,g,g + g') ~ g^A^ . 


(28) 
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This is by a factor {qy/A)'^ smaller than the unsym¬ 
metrized three-point loop, that is, the degree of UV di¬ 
vergence has been reduced by two upon symmetrization. 
The sum over the two inequivalent Aslamasov-Larkin dia¬ 
grams is thus ultraviolet-finite at least for z < 4, as shown 
previously for z = 3 by Metlitski and Sachdev [5]. 


V. ASYMPTOTIC PROPERTIES OF Y-POINT 
LOOP 


is an antisymmetric polynomial of the form 


A^ijlfijl,u (^o) 


Pix Pjx 

Pix 

Pvx 



Piy Pjy 

Ply 

Pvy 



P% P% 

Pfy 

Pvy 




1 1 

1 

1 




Pix 

Pjx 

Pi 

r 

Pux 


Piy 

Pjy 

Ply 

Pvy 

z 

{ho} 

{kjo} {ho} 

{kvo} 


1 

1 

1 


1 


(32) 


We have seen that the three-point loop is cut off by 
a kinematic constraint at large momenta, if one of the 
external momenta stays fixed, and cancellations occur 
upon symmetrization. In this section we investigate to 
what extent these properties can be generalized to N- 
point loops. To this end, we first rewrite In{pi, ■ ■ ■ ,Pn) 
in a form which is more convenient for an asymptotic 
analysis. 


with kio = ko - Pio- Hence, Dijifiji^^{ko) is totally anti¬ 
symmetric in all four indices. Since Diji is also antisym¬ 
metric, it follows that is invariant under permuta¬ 
tions of i,j, 1. 

For loops constructed with bare propagators Go, the 
quantities Jiji and fiji^u are independent of fco- The 
frequency integration in Eq. (30) can then be carried out 
and one recovers the reduction formula for /^(pi,... ,pAr) 
in the form Eq. (29). 


A. Reduction formula for Y-point loop 


B. Kinematic constraint for Y-point loop 


Eeldman et al. [17] derived a reduction formula ex¬ 
pressing the Y-point loop with bare propagators and a 
quadratic dispersion relation as a linear combination of 
three-point loops, 


I%{pi,...,Pn) = f,Y ^ ( 29 ) 




where fiji,v is a rational function oi pi,pj,pi,pi, [18]. 

For the Y-point loop with non-Fermi liquid propagators 
of the form Eq. (5), a similar formula can be derived. In 
Appendix B we show that In{pi, ■ ■ ■ ,Pn) as given by 
Eq. (11) can be expressed as 


We now analyze the asymptotic behavior of an Y-point 
loop traversed by a single large momentum. More pre¬ 
cisely, we assume that two of the external momenta, say 
q' and q" k, —q' are large and almost antiparallel, while 
all the other momenta are kept finite, that is, relatively 
small. This implies that the momenta Pi fall in two groups, 
where momenta within a group are close together, while 
the distance between momenta in different groups is large 
(see Fig. 3). Choosing pi — 0, the momenta in the group 
containing pi are all close to zero, that is, relatively small. 

Using the generalized reduction formula (30), we can 
show that the frequency integrations are effectively re¬ 
stricted to small intervals, although the differences Pio—pjo 
may be large. To this end, let us fix some i^j, I and write 


E / (r®( 


PiO — PjO 

Piy — Pjy 


n 4/.i'(^o) Jiji{ko) +eye, (30) 


where 




Dijl{Piy ~ Pjy) 


(31) 



with Jiji{ko) and Diji from Eqs. (12) and (13), respec¬ 
tively. “eye” denotes cyclic permutations of i,j, Z. Ex¬ 
pressing Jiji{ko) and Diji explicitly in terms of momenta 
and frequencies, one finds that the product Dijifiji^,^{ko) 


Figure 3. Y-point loop with two large momenta q' and q" ~ 
—q'. Momenta pi on the bold lines are large, those on thin 
lines (relatively) small. 
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out the cyclic sum: 


+ 

+ 



PiO PjO 

Piy — Pjy 

PjO — Pm 
Pjy ~ Ply 
PlO — PiO 
Ply Piy 




(33) 


where 5 (^ 0 ) = I\u^i,j,i fiji,ui^o) Jzji{ko)- Note that the 


latter function is invariant under cyclic permutations of 

If Pi, pj, Pi are in the same group of small or large 
momenta, all frequency integrations in (33) are obviously 
limited to small intervals. Now assume that they are in 
two distinct groups, say Pi is large and the others are 
small. The argument of the first and third step function 
in (33) is then dominated by pi, leaving Q{pio/piy) in 
both cases, and the whole expression can be simplified to 


'PjO 


dko 

47r 


0 



dko ^f PjO -Pio 
Pjy ~ Ply 


Jv,n 47r V 


(34) 

The frequency integration is thus effectively restricted 
to the small interval between pjo and piQ. The large 
contributions from integrations from pjQ to PiQ and from 
Pio to pio cancel in the cyclic sum. 

In naive power-counting for the 7V-point loop, the fre¬ 
quency integration yields a factor A^. In the specific 
limit discussed above, the ultraviolet asymptotics is thus 
reduced by that factor. For the above argument it was im¬ 
portant that both Qq and q'y are large. The ^-component 
did not matter. If only q'g is assumed to be large, the 
/cQ-integration does extend over a large interval, yielding a 
factor A^, but only as long as q'y is of the order of the other 
y-components. One thus obtains an effective restriction of 
q'y as discussed already for the three-point loop in Sec. IV. 

In summary, we have shown that a fermion loop tra¬ 
versed by a single large momentum is suppressed by cancel¬ 
lations in the cyclic sum in Eq. (30). No such cancellation 
occurs if three or more momenta qi are large. In that case 
three momenta pi , pj , pi can be far apart from each other 
so that no systematic cancellation occurs in the cyclic 
sum. 


C. Symmetrized N-point loop with one or two 
small momenta 

We now show that systematic cancellations occur in 
the symmetrized N-point loop in the limit of one or two 
vanishing momenta. This generally reduces the degree of 
divergence of Feynman diagrams with bosonic external 
legs coupled to a fermion loop. 

We first show that the symmetrized iV-point loop van¬ 
ishes, if one of the external momenta vanishes. The sym¬ 
metrized iV-point loop is given by a sum over all permu¬ 
tations of external momenta. For a vanishing external 


momentum two fermion lines in the loop carry the same 
internal momentum. By permutations the vertex with the 
vanishing momentum is cycled around the loop, yielding 
the sum 

/* rlh r 1 N—1 

j ^ j ( 35 ) 

Since the denominator of G{k — Pi) is linear in 
kx, the integrand can be written as a /Ca;-derivative, 
— g|- G{k—pi). Performing the fca;-integration one 

thus finds that the symmetrized iV-point loop vanishes 
if one leg has a vanishing momentum. By dimensional 
analysis, the UV scaling of the symmetrized loop with 
one fixed external momentum is thus reduced by a factor 

A-i. 

Using analyticity and the invariance under qi 1 —> —qi 
one can conclude that symmetrized loops vanish even 
quadratically, if a vanishing momentum q enters and 
leaves the same loop at two distinct vertices, provided 
that the other momenta remain finite. Due to momentum 
conservation this is possible only for A^ > 4. Hence, the 
UV scaling of symmetrized loops for N — 2 large momenta 
and two fixed external momenta q and —q is reduced by 
a factor A“^. 


VI. IMPROVED POWER-COUNTING 

We now use the above results on the asymptotic behav¬ 
ior of fermion loops to obtain improved power-counting 
estimates for the ultraviolet behavior of several classes 
of Feynman diagrams. We consider contributions to the 
boson and fermion self-energies, the fermion-boson vertex, 
and the three-boson vertex. 


A. Boson self-energy 

Perturbative contributions to the boson self-energy 
T,b{q) with L loop-integrations contain L — 1 boson prop¬ 
agators and 2L fermion propagators. The boson prop¬ 
agators decay as A^“^ for large momenta, the fermion 
propagators as A~^. Each loop integration contributes a 
power Hence, according to naive power-counting, L- 

loop contributions to the boson self-energy may diverge as 
A^~^. In particular, for z = 3 a quadratic UV divergence 
seems possible for any loop order L. 

However, systematic cancellations of UV divergences oc¬ 
cur between distinct contributions (represented by distinct 
Feynman diagrams) at any given loop order. Feynman 
diagrams contributing to the boson self-energy ^tiq) con¬ 
tain fermion loops but no open fermion lines. Hence, 
the external boson legs with ingoing and outgoing mo¬ 
mentum q connect directly to fermion loops. The sum 
over all diagrams of a certain order can be written in 
terms of symmetrized fermion loops. In Sec. V.C we have 
shown that in symmetrized fermion loops with a fixed 
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Figure 4. A class of four-loop diagrams with fermion loops 
connected to three boson propagators with possibly large 
momenta. 

(relatively small) momentum q, the UV divergences are 
suppressed by a factor upon symmetrization, and 
in symmetrized loops with two hxed momenta q and —q 
by a factor A“^. In Feynman diagrams contributing to 
the boson self-energy the external momenta q and —q are 
either attached to two distinct loops or to just one loop. 
In both cases there is a power-counting gain of the order 
A~^. Hence, the sum over all L-loop contributions to the 
boson self-energy behaves as for large A. For z < 3 
this means that the sum is UV-convergent. For z = 3, 
quadratic and linear UV divergences must cancel. This 
cancellation is also imposed by a Ward identity following 
from current conservation [5]. However, logarithmic UV 
divergences are not excluded for z = 3. 

At three-loop order, the only candidates for a logarith¬ 
mic UV divergence for z = 3 are the Aslamasov-Larkin 
diagrams (see Fig. 2). However, as discussed already in 
Sec. IV, these diagrams are restricted by a kinematic 
constraint arising from cancellations in the representation 
of the three-point fermion loop as a cyclic sum, so that 
the total contribution is UV-finite for any z < 4. 

As shown in Sec. V.B, a kinematic constraint from a 
cancellation of contributions to the V-point fermion loops 
in the UV limit occurs only if only one large momentum 
traverses the loop. Contributions to the boson self-energy 
where this condition is not satisfied occur at four-loop 
order. An example is shown in Fig. 4. In Ref. 10, those 
contributions (with symmetrized fermion loops) were com¬ 
puted explicitly for z = 3 in the static limit qo = 0, and 
were indeed found to diverge logarithmically. Surprisingly, 
the divergence turned out to be of the order log®(A/|gj^|) 
instead of the expected simple logarithm. Such a diver¬ 
gence is not renormalizable. The implications of that 
singularity are not clear at the moment. 

It is remarkable that the one-loop result receives quali¬ 


tative corrections only at four-loop order, while the dy¬ 
namical exponent z = 3 remains unchanged at two- and 
three loop order. Such a situation is peculiar but not 
unprecedented in quantum field theory. For example, in 
non-linear tr-models describing the critical behavior of 
the Anderson localization transition, only the four-loop 
contributions lift a degeneracy between distinct symme¬ 
try groups and correct the one-loop result for critical 
localization length and conductivity exponents [19]. 


B. Fermion self-energy 

Perturbative contributions to the fermion self-energy 
E/(fc) with L loops contain L boson propagators and 
2L — 1 fermion propagators. Hence, L loop contributions 
to the fermion self-energy diverge as A^ in the ultraviolet 
limit for any L and z. This divergent contribution is 
however momentum and frequency independent and can 
be absorbed by a shift of the Fermi surface. Physically 
more interesting is the momentum and frequency depen¬ 
dence of the self-energy. In particular, the derivative 
f /9/c+ with k+ = kx + ky determines the anomalous 
dimension of the fermion fields [5]. Since kx and ky scale 
as A^, that derivative scales as A° according to the above 
power-counting. Hence, logarithmic UV divergences are 
expected, and indeed occur already at three-loop order 
(see Fig. 5), as discovered for z = 3 by Metlitski and 
Sachdev [5]. The three-point fermion loop in these three- 
loop contributions is not reduced by any of the cancella¬ 
tions we discussed above, since all three boson momenta 
can be large. 



Figure 5. Logarithmically divergent three-loop contributions 
to the fermion self-energy. 


C. Fermion-boson vertex 

We now discuss the fermion-boson vertex r(fc, q) with 
one boson and two fermion legs. Perturbative L-loop 
vertex corrections contain L boson propagators and 2L 
fermion propagators. Hence, F(fc, g) scales as A° in the 
UV limit, that is, logarithmic divergences are expected. 
The vertex corrections can be grouped in two classes. In 
the first case, the external boson leg couples to an open 
fermion line, and in the second to a fermion loop. In Fig. 6 
we provide a three-loop example for each case. In the 
first case we do not expect any cancellations, but in the 























Figure 6. Two types of vertex corrections with the external 
boson leg linked to an open fermion line (left) and to a fermion 
loop (right). 
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Figure 7. Four-loop diagram contributing to the three-boson 
vertex. 


second we have a fermion loop with a fixed boson momen¬ 
tum q. Hence, summing all diagrams corresponding to 
the symmetrized fermion loop, the UV contributions are 
suppressed by a factor so that the sum is guaranteed 
to be finite. 

The fermion-boson vertex is related to the fermion self¬ 
energy by the usual Ward-identity following from charge 
conservation [5]. Since the fermion self-energy exhibits a 
logarithmic divergence already at three-loop order, the 
fermion-boson vertex has to diverge, too. The above 
argument shows that the divergent contributions to the 
vertex come exclusively from those diagrams where the 
boson leg couples to an open fermion line. 

D. Three-boson vertex 

As a final example, we consider the three-boson vertex. 
Its lowest order (one loop) contribution is simply the 
symmetrized three-point fermion loop Hg^™, which does 
not exhibit any UV divergence. L-loop corrections to the 
three-boson vertex contain L — 1 boson propagators and 
2L + 1 fermion propagators. Naive power-counting thus 
yields a UV scaling of the form A^“^. For z = 3 this 
is marginal. However, each of the external boson legs is 
linked directly to a fermion loop (see, for example. Fig. 7). 
Upon loop symmetrization possible UV divergences are 
thus suppressed. Without investigating the exact degree 
of suppression we can say that in the most interesting 
physical case z = 3 all corrections will be UV finite. The 
UV convergence of Wloop vertices with > 3 is also 
improved by cancellations in the symmetrized fermion 
loops, but they are already UV convergent even within 
naive power-counting for z < 2N — 3. 

VII. CONCLUSION 

We have analyzed general properties of the perturba¬ 
tive loop expansion for two-dimensional quantum critical 
metals with singular forward scattering. Important cases 
falling in this class of systems are metals at an Ising ne¬ 
matic QCP and metals coupled to a U{1) gauge field. 


Our analysis is based on the effective field theory for such 
systems [5,9], extended to arbitrary dynamical exponents 
z > 2 as proposed earlier by Nayak and Wilczek [11]. 

We have derived asymptotic properties of fermion loops 
appearing as subdiagrams of Feynman diagrams contribut¬ 
ing to the perturbation series. Substantial cancellations 
occur when summing over contributions corresponding 
to permutations of vertices at fermion loops with a small 
external momentum q or two small momenta q and —q 
(while all the other momenta become large). Further¬ 
more, fermion loops are suppressed below the naive power¬ 
counting estimate when they are traversed by only one 
large bosonic momentum. 

Using these properties we have shown in a number of 
specific examples how one can perform a sharp ultraviolet 
power-counting that takes all possible cancellations into 
account. Most importantly, for the boson self-energy there 
is always a gain of order A^. As a consequence, perturba¬ 
tive contributions to that quantity are UV convergent for 
z < 3. Hence, perturbative self-energy corrections to the 
boson propagator are finite at all loop orders for z < 3. 
In the static limit qo = 0, the boson self-energy is then 
fixed by dimensional analysis to be proportional to qy~^ 
with a finite prefactor. A bare z < 3 is thus not renor¬ 
malized by fluctuations. On the other hand, logarithmic 
divergences are possible in the important physical case 
with a dynamical exponent z = 3, and indeed do occur at 
four-loop order [10]. For z > 3 there are UV divergences 
proportional to A*“^. 

Obviously, z = 3 is a critical value at the boundary 
between the qualitatively different cases z < 3 and z > 3. 
Interestingly, the same critical value z = 3 separates two 
qualitatively different behaviors of the compressibility, if 
the chemical potential /i couples to the bosonic mass ap¬ 
pearing away from the QCP. For z < 3 the compressibility 
is hnite, while for z > 3 it diverges upon approaching the 
QCP [5]. Our result that a bare z < 3 is not renormalized 
by fluctuations at any loop order thus implies that the 
compressibility remains finite in that case. 

There are no systematic cancellations reducing the 
UV divergences of the fermion self-energy. Hence, the 
fermionic field renormalization obtained from a first mo¬ 
mentum derivative of the self-energy diverges logarith- 
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mically for any z. In contrast to the behavior of the 
boson self-energy, the degree of divergence of the fermion 
self-energy is independent of z. Logarithmic divergences 
obtained in perturbation theory can be summed up to 
yield anomalous power-law scaling, as demonstrated in 
the three-loop calculation for z = 3 by Metlitski and 
Sachdev [5]. 

At this point, the perturbative singularity structure of 
the above class of critical metals seems to be clarified for 
z < 3. The one-loop fixed point is only modified by an 
anomalous scaling of the fermion fields. However, the final 
fate of the theory in the special but physically relevant 
case z = 3 remains open. Organizing the loop expansion 
in terms of symmetrized fermion loops and exploiting the 
cancellations described above should be useful for any 
future work on fluctuation corrections. In particular, a 
numerical computation of fluctuation corrections is greatly 
facilitated by incorporating the cancellations directly via 
symmetrized loops. 
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Appendix A: Integration of loops by residues 


Here we derive Eq. (11) for the A^-point fermion loop 



(Al) 

where G'(fc) = \kx + — i{fco}] ^ with {ko} = fco/|A:o|“. 

The fca;-integration can be performed by closing the inte¬ 
gration contour in the upper complex half-plane and using 
the residue theorem. There are N simple poles situated 
at kx = i{ko - Pm} + ptx - {ky - piy^ with * = 1,..., A. 
Summing the contributions from all poles in the upper 
half plane (fco > pm) yields 


IM{Pl,■.■,PN) = ^ ^ J_ 

xn 


dky 

27T 


N 


^ 0 (/co - Vm) 
1 


i{ko - Pm} - i{ko - Pjo} + - Pjx + ‘^iPty - Pjy)ky - pf + p] 


(A2) 


'L-u rtyjj -L-u x-juj • nu. jrju. • -\r-Ly Jr j y / ■ ~ y r'ly ' t' jy 

Now the A:y-integration can be performed analogously. For each i there are N — 1 simple poles at 


ky = 


i{ko - Pjo} - i{ko - Pm} - Pix + Pjx + pfy - p]y 
‘^{Piy ~ Pjy) 


with j 7 ^ i. These poles are situated in the upper complex half plane if (and only if) > 0. Closing the integration 

Piy Pjy 

contour in the upper complex half plane thus yields 


Defining 


In{pi,---,Pn) = - 


dko 

27T 


Y^'^Qiko -Pm)Q(^ 

i j^i 


PiO — PjO \ _ 

Piy — Pjy^ ‘^{Piy ~ Pjy) 


X It [^{^0 - Pio} - i{ko - Plo} + Pix - Plx - Pfy + Pfy 

+ ——^(*{^0 - Pio } - i{ko - Pm} - Pix + Pjx + P% - Pm) 
Piy -Pjy 


Diji = -det 


Pix Pjx Pix Plx 
Piy — Pjy Piy — Ply 


= PixiPly - Pjy) + eye , 


Fiji — {Pjy Piy){Ply Pjy)iPiy Ply) ^ 
^ijiiko) = {ko -PM}iPiy -Pjy) + cyc, 

Eq. (A3) can be written as 


(A3) 


(A4) 

(A5) 

(A6) 


In^Pi 


,Pn) = - 


i j^i 


PiO — PjO 
Piy — Pjy 


(Piy Pjy) 


N-3 


It [Diji + Fiji 

l9^ij 


-iiliji{ko)] (A7) 


Defining Jijiiko) = [Diji + Fiji + iU,iji{ko)] and using the antisymmetry Jiji{ko) = —Jjii{ko), one obtains 
In{pi, . ■ ■ iPn) as a sum over indices i and j restricted to i < j in the form Eq. (11). 
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Appendix B: Derivation of reduction formula 


In this Appendix we derive the expression Eq. (30) for the A-point loop from Eq. (11). The crucial step is the 
identity 


where fiji,v{ko) is defined as in Eq. (31), that is, 


fijl,viko) — 


{piy Pjy) Jijliko) — fijl,iy{ko) Jijl{ko) . 


(Bl) 




^ijl^ijuiko) ^ijvJijl (^o) 
DijliPiy ~ Pjy) 


(B2) 


This can be verified by first considering the cases A = 4 and A = 5, and then proceeding by induction. For A = 4, 
Eq. (Bl) reads (suppressing the fco-dependence) 


[Piy Pjy)JijliJijli — J 


(B3) 


which follows from 


JijhJijh ~ '^ijh^ijhi^ 


Dijli Jijli 


+ ll 1 - 2 ^ — 


D, 


ijh 




rziJijh W I 2 ■ 


' DijhJijh ^ijhJijh 
For A = 5, one finds from repeated application of Eq. (B3) 

(Piy — Pjy) JijhJijl^Jijh — /yii ,Z2 h I3 ^■^ijhji-l^ijh.h 

.^q7i ,/2 I 3 *^^ 772 li'^q72,^3T 


(B4) 


where in the last step we have used = 44 ./i 4 bl 2 ’ ^e verified by explicit 

calculation. It is now clear that Eq. (Bl) follows from (A^ — 3)-fold iterated application of Eq. (B3). 

Inserting Eq. (Bl) into Eq. (11) yields 


n-l n -1 


n -1 p -1 


/»(P..p») = I:E /""fse)™^) 




47r \p^y-Pjy^ 


H hji,v{ko) Jiji{ko). 


(B5) 


Using the invariance of Jiji{ko) and fiji^viko) under cyclic permutations of i,j, I one obtains In{pi, ■ ■ ■ ,Pn) in the 
cyclic form Eq. (30). 
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